Intrinsic Spin-Orbit Interaction in Graphene 
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In graphene, we report the first theoretical demonstration of how the intrinsic spin orbit interaction 
can be deduced from the theory and how it can be controlled by tuning a uniform magnetic field, 
and/or by changing the strength of a long range Coulomb like impurity (adatom), as well as gap 
parameter. In the impurity context, we find that intrinsic spin-orbit interaction energy may be 
enhanced by increasing the strength of magnetic field and/or by decreasing the band gap mass 
term. Additionally, it may be strongly enhanced by increasing the impurity strength. Furthermore, 
from the proposal of Kane and Mele [Phys. Rev. Lett. 95, 226801 (2005)], it was discussed that 
the pristine graphene has a quantized spin Hall effect regime where the Rashba type spin orbit 
interaction term is smaller than that of intrinsic one. Our analysis suggest the nonexistence of such 
a regime in the ground state of fiat graphene. 

PACS numbers: 73.22.Pr,71.70.Ej,73.43.Cd,72.80.Vp 



The spin-orbit effect is purely kinematic effect, and 
arises due to the interaction of the orbital and spin inter- 
nal degrees of freedom. Recently, there has been growing 
interest in the study of this effect in graphene, due the 
fact that the progress in manipulating graphene^, and 
graphene based nanostructures, particularly graphene 
quantum dots^^— has opened new perspectives in usage of 
spin qubits^ for quantum computing and quantum infor- 
mation purposes^. Graphene based nanostructures seem 
to be proper candidate for realizing these kind of pro- 
cesses through controlling one or two basic spin inter- 
actions such as intrinsic or Rashba type spin-orbit in- 
teractions (SOIs)''"^^. Since the former is supposed to 
be small, investigations are focused on the latter one 
which assigns a finite mass for the graphene carriers 
due to the interaction with substrate, and thus lead- 
ing to a finite gap. In fact, within the Spin Quantum 
Hall Effect (SQHE) context, investigations were initi- 
ated by Kane and Mele's worki. They proposed that 
the ground state of graphene exhibits a SQHE. They 
showed that SQHE regime exists provided that Rashba 
type SOI term, AgQ ~ 0.05 meV, smaller than that 
found for intrinsic SOI, AgQ ^ 0.103 meV. Later, on 
the basis of microscopic considerations, AgQ (AgQ ) is es- 
timated to be much smaller (larger) than that of found 
by Kane and Mele: For instance, Hernando et aU^ found 
A|q - 0.86 /xeV, Min et al2. found A^gQ ~ 0.57 /leV, 
AgQ ^ 0.011 meV. Nevertheless, it is concluded that 
SQHE regime can occur only below a temperature AgQ ~ 
0.01 K (10""^ meV). In the impurity context, very re- 
cently Castro and Netoi2, showed that impurity induced 
distortion in flat graphene lead to a significant enhance- 
ment of SOI about AgQ^ ^ 1—7 meV. While our analysis 
confirms this estimate within the pure impurity coverage, 
but it reveals that, in the absence of impurities, i.e., in 
the case of flat graphene, there is no experimentally ac- 
cessible temperature where SQHE regime can occur. 

Since intrinsic SOI is a natural property of Dirac's 
equation, i.e., it inherently involves the spin degrees of 
freedom, its Graphene's analog contains pseudospin or 
just spin-orbit interaction automatically. By graphene's 



spin analog, we refer two equivalent points in the Bril- 
louin zone, i.e., K and K valleys as pseudospin or just 
gpij^20_ this letter, we flrst propose a simple model 
to address the question of how intrinsic SOI energy in 
graphene can be obtained, and then can be controlled by 
the strength of applied magnetic fleld, and/or impurity 
coverage as well as gap parameter. 

In graphene, near the Dirac point, the electronic states 
in the presence of electromagnetic potential = 
(Ao,A) are described by the effective low-energy Dirac 
equation, HvP = E'^ where 



"H = vpOL • ( p A 



eAo 



(1) 



is the Dirac Hamiltonian. Here, A and Aq are the vec- 
toral and scalar parts of the four vector potential A^, re- 
spectively, and they will be chosen as A = B(— y, x,0)/2 
in the Goulomb gauge, and Ze/er for charged Coulomb 
impurity, respectively, oti and f3i are Dirac matrices, 
and vp ~ lO^m/s is the Fermi velocity. In the ab- 
sence of external fields, the energy eigenvalues of T-l yield 
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E = ±\/ Mq + /cq which determines linear dispersion 
ko (~fco) for the conduction (valence) in the absence of 
a gap term Mq — gMq = 2 A/3. We have used abbrevi- 
ations Mq — A4vp/fi = 2A/3a together with definitions 
fivp = 3aJo/2, where a and Jq are the carbon-carbon dis- 
tance (1.42A) and transfer integral (~ 2.7eV) between 
them, respectively. Thus, E = Ea/Tivp and k = ka are 
dimensionless energy and wave vector. Throughout the 
letter, we restrict ourselves in a single valley ( K ) and in 
a single band (conduction). 

By successively applying Fouldy-Wouthuysen unitary 
transformation on wave functions and operators of Dirac 
equation 'H'^ = i?^', one passes the two component equa- 
tion at any desired order of electromagnetic interaction 
strength, and thus the small and large components of ^ 
are completely decoupled^'^. Hence, by just restricting 
ourselves to the positive energy solutions, together with 
the corresponding upper components of ^, expansion of 
% in power of the strength of electromagnetic interaction. 
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as it should be, yields the well-known two component and j = ruj + 1/2 is the eigenvalue of the conserved to- 

Lz + Sz. The values of the 
if rrij > 0, and 



Pauh equation, = (Mvj. + H„.i,)$ where 



c 



eh 
2mc 



a- B-eAo (2) 



whose eigenvalue equation Hn.i.^ — yields well- 

known nonrelativistic results, but in the two components 
formalism. Continuing the expansion to second order, 
one can obtain three more terms that we represent them 
as A, thus, -H^w^ ^ _^ ^. ^ ^-^^ rpj^ggg 

relativistic dependence of kinetic energy, SOI energy, and 
the Darwin terms, respectively. The last term is different 
from zero at where there are charges creating the field^^. 
Furthermore, it is ascribed to ZitterbewegungSi., and may 
be dominated at higher magnetic fields. Thus, just by 
taking the difference of energy spectra of H and 7{^^ 
given by Eq. ([T]) and Eq. ([2]), respectively, we are left with 
SOI energy. In other words, it may be defined as Aii^so = 
—Mvp + E — En.r. , where E is the relativistic energy 
spectrum. In dimensionless form, it can be rewritten as 
AEso = -Mo + E-En.r.. 

We first discuss the influence of magnetic field on the 
graphene electronic energy spectrum in the presence of a 
single Coulomb impurity. In this case, the Dirac Hamil- 
tonian is by 



where 



Ho 



H = 'Hq — fivp-a- ■ A, 

c 



ihvF (cr • V) - {Ze^/er) 



(3) 



is the exactly solvable unperturbed part^*^ whose en- 
ergy eigenvalues are easily found to be as e„ = 



M 



—2 



-1/2 



1 + Z / (ti + 7) together with the corre- 

sponding eigenfunctions in terms of Laguerre polynomi- 
als 



]_ ( (r) e'O-i) 

V I Gnj (r) e'(^) 



(4) 



where 



with 



J^] I = (-l)"iV,y (^,Mo) ± e„ 



-Ar 



(2Ar) 



7-1/2 



L^,^ (2Ar) ± (2Ar) 



{Z, Mo) 



r{n+l)X^ [rUj + (ZMn/X)] 



1/2 



r (n 27 -f 1) ZM§ 



is the normalization constant, and 



C2 



(n-f 27)/[mj-f (ZMo/A)], A = 



0-4, 7 



-Z , and Ll^ are Laguerre polynomials. Here, 
Ze'^/ehvp is the dimensionless coupling constant 



tal angular momentum J; 
quantum number n are n = 0,1,2,. 
71 = 1, 2, ... if nij < 0. It can be easily check that eo — 



Mo 



1 



2Z) 



1/2 



becomes zero at Zcr — 1/2. In the 

framework of perturbation theory, by using Eq. Q, one 
obtains first-order shift in energy eigenvalues of Eq. ^ 
as E = In + Ae„ , where 



1 



4Z 



B^Mo~el{2n + ^)j + 



ZMo 



'Mo-n 
(5) 

is the Zeeman term due to the magnetic field. In Eq. ([5|) , 

B = B/Bq is the dimensionless magnetic field with 

Bo = fic/ea^. The result given in Eq. ([5]) is valid for 
2 2 

B <^ Z Mq/2. In the same units, the energy eigenval- 
ues En together with its two-dimensional non-relativistic 
equivalent which can easily be obtained^^ as 



E„ 



1 



B Mo 



(m -|- s) 



can be used to calculate SOI energy as functions of 
strengths of magnetic field, gap A parameter as well as 
of impurity strength Z. 

First, in the absence of magnetic field, to show impu- 
rity strength dependency of the SOI energy alone in the 
graphene, we plotted the lowest lying states in fixed gap 
values as a function of Z in FIG. [TJa). For comparison, 
in the inset, we have also included a plot showing the gap 
parameter dependency of SOI energy for the same energy 
levels but at fixed impurity strengths. It is clear from the 
figure that, in comparison with SOI energies for higher 
orbital, SOI energy for the ground-state is higher, and 
all they are getting more pronounced with increase in Z 
and/or A. The former follows from the fact that, higher 
states signify larger radius so that the corresponding im- 
purity binding energy, and hence the SOI energy becomes 
smaller. The latter can easily be approved by analyzing 
the coupling constant and mass dependent form of SOI 
energy for two-dimensional relativistic Hydrogen atom2£. 
As a result, not only our these findings support the en- 
hancement of SOI by impurity covarag&i^ , but they also 
show how it can be controlled by means of the relevant 
parameters such as impurity strength and gap parameter. 

In the presence of uniform magnetic field, as is seen 
from Eq. ([S]), the orbital degeneracies are lifted. In par- 
ticular, it splits the first excited doubly degenerate j lev- 
els, and thus leads to Zeeman splitting. By increasing the 
magnetic field, as expected, the Zeeman splitting spac- 
ing broadens gradually with B, (FIG.HJb)). To show the 
combined effects of magnetic field on SOI together with 
the effect of SOI on Zeeman splitting, we plot the Zee- 
man splitting for the first excited state in the relativistic 
spectrum as function of magnetic field in FIG. [Ifb). We 
also include a corresponding curve for the non-relativistic 
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FIG. 1. (Color online) (a) SOI energies as function of Z at 
gap parameter A = 0.1 (straight lines) and A = 0.01 (dashed 
lines) for the three lowest-lying energy levels, namely, from 
top to bottom, [n,j] = [0, 1/2], [1, -1/2] = [l,+l/2] and 
[0,-3/2] = [0,3/2], respectively. Inset shows its dependence 
on A at Z = 0.1 (dashed lines) and Z — 0.4 (straight lines), 
respectively, for the same levels, (b) Relativistic (bold black 
line) and non-relativistic (dashed black line) Zeeman splitting 
energies in [1,-1/2] , [1,-1-1/2] levels, together with their dif- 
ference (red straight line) as function of magnetic field at fixed 
value Z — 0.3. For this value of Z , we also plot individual SOI 
for each level, i.e., [1,-1/2] (green), [1, +1/2] (Blue) as func- 
tion of magnetic field. Inset: relativistic (dashed lines) and 
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energy spectrum as a reference. Hence, by just taking 
the difference between these two, i.e., between the rela- 
tivistic and the non-relativistic energies Zeeman energies, 
we obtain the contribution of SOI energy to the Zeeman 
splitting energy. Its dependence on B, magnetic field de- 
pendence of SOI for j = 1/2 and j — —1/2 levels are 
individually displayed in FIG. lUb) for comparison. We 
can see that its influence becomes significant for decreas- 
ing B, since the Zeeman splitting linearly depend on B 
and j. Therefore, the difference between these two curves 
is simply the SOI contribution to the Zeeman energy. 
To learn the response of SOI to the applied magnetic 
field, we have also plotted the individual SOI energies 
for each level. Obviously, we claim that, when the asso- 
ciated Zeeman splitting is resolved in an experiment, the 
related SOI contribution, and therefore SOI energies for 
each level can be resolved. As an example, ai B — IT, 
our calculation describes Zeeman energy up to 0.70 meV 
which contains O.lSmeV due to SOI. One can also jus- 
tify this value by just taking the difference between the 
SOI energies of the related levels which are 0.20 meV and 

0. 05 meV , respectively. As for the inset of the figure, we 
compare A dependency of the SOI energies. Eventually, 
we give a measure of how amount of Zeeman energy is 
influenced by SOI. The above values may be enhanced 
by just increasing Z and/or A decreasing . 

A similar procedure developed above for binding in 
the hydrogenic sense can also be extended to the binding 
in oscillatory sense, in the absence of any impurity. In 
pristine graphene, massless Dirac Fermions having linear 
dispersion obey Dirac- Weyl equation, and they have no 
nonrelativistic analog. But, we may overcome this diffi- 
culty by just using the spectrum of Dirac equation, and 
then take its limit when A — >■ 0. To do this, it is enough 
to recall the well-known graphene Landau levels (LLs), 

1. e., 

E^\JmI + 2n/4 

with n = i^+{\m\ +m + s + 1) /2. Here, is = 1/Vb is 
the magnetic confinement length. It can easily be shown 
that its non-relativistic counterpart has energy spectrum 
with 

— 1 / ^ \m\+m+l \ g*5 

^n.r. = 2 ^ o + 2 

wherein the first term is the well-known non-relativistic 
LLs, the second one is due to Zeeman effect. Here, we will 
set the Lande factor to 2. It is easy to check analytically 
that SOI energy is zero for n = {v = 0,m = 0, s = — 1) 
level even without taking A — > limit (massless limit). It 
is also easy to see that higher orbital are sensitive to de- 
crease in A, i.e., A — limit (gapless limit). As is done 
in FIG. [21 it is also instructive to compare this picture 
with those obtained from the free particle spectrum, i.e., 

A^so -'Mq + \Jl^l + fco - (^A:o/2Mo) which is also 
zero at fco ~^ limit. Again, even in the presence of mass 
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FIG. 2. (Color online) For the four lowest eigenstates of 
the graphene, SOI energies as function of gap parameter at 
B = 1 T. The upper scale shows the gap parameter in 
units of meV. Inset: SOI for free particle spectrum, i.e., 

^Eso = —Mo + \J All + fco — ^feo/2Mo^ as function of gap 

parameter. From top to bottom, red, green and black straight 
lines corresponds fco — 0.1, 0.01, and 0.001, respectively. 



term, it gives zero. Whereas it becomes more pronounced 
as A — > for finite ko. The appearance of the lack of SOI 
in the electron- hole degeneracy point, i.e., in both n = 
and fco = cases, shows that this point has purely non- 
relativistic character. Furthermore, this proves that the 
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ground-state of the flat graphene does not exhibit SQHE 
regime. Additionally, beyond this point we see that SOI 
effects are much more pronounced in A — > limit than 
those found for the impurity coverage in the same limit. 
These findings are exactly compatible with the analyti- 
cal predictions of the relativistic quantum mechanics on 
the general form for SOI energy in general. All these 
show that our predictions work well for different physical 
regimes. 

We have three main conclusions. First, we demon- 
strate theoretically that, without sophisticated many- 
body calculations, intrinsic SOI energy in the electron- 
hole degeneracy point is zero in flat graphene. Thus, we 
conclude that, for flat graphene, it is not possible to ob- 
serve SQHE. Second, within the impurity coverage, we 
investigate the influence of SOI on graphene energy lev- 
els. Finally, in the same context, we investigate the Zee- 
man splitting in doped graphene as well as the influence 
of SOI on this splitting . We demonstrate that it is possi- 
ble to enhance intrinsic SOI by increasing the strength of 
applied magnetic field, and how it can be resolved from 
the Zeeman effect. Therefore, magnetic field not only 
can be used to tune the Zeeman effect, but also allow us 
to estimate the magnitude of the SOI energy in dopped 
graphene. As for the Zitterbewegung, its effect should 
be taken account for higher levels and higher magnetic 
field values. Its signature is the oscillatory behavior of 
the velocity and for its direct observation stronger mag- 
netic fields would be required^^. It can also be resolved 
in the same procedure described here within this con- 
text, i.e., but for higher levels where SOI is negligible, 
and for higher magnetic fields. We hope that this study 
will not just reveal the intrinsic SOI in graphene, but 
also it will illuminate the Berry phase in graphene, since 
very recently it is shown that, by constructing the intrin- 
sic SOI operator from the non-Abelian Berry connection, 
intrinsic SOI originates from Berry phase terms^^. 
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